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I. INTRODUCTION 



The phase-space formulation of quantum mechanics 
has a long history. In 1932 Wigner yj] introduced his 
famous function which has found numerous applications 
in many areas of physics and electronics. In 1949 Moyal 
H discovered that the Weyl correspondence rule || can 
be inverted by the Wigner transform from an operator 
on the Hilbert space to a function on the phase space. 
As a result, the quantum expectation value of an op- 
erator can be represented by the statistical-like average 
of the corresponding phase-space function with the sta- 
tistical density given by the Wigner function associated 
with the density matrix of the quantum state. In this 
way quantum mechanics can be formally represented as 
a statistical theory on classical phase space. It should 
be emphasized that this phase-space formalism does not 
replace quantum mechanics by a classical or semiclassical 
theory. In fact, the phase-space formulation of quantum 
mechanics (also known as the Moyal quantization) is in 
principle equivalent to conventional formulations due to 
Heisenberg, Schrodinger, and Feynman. However, the 
formal resemblance of quantum mechanics in the Moyal 
formulation to classical statistical mechanics can yield 
deeper understanding of differences between the quan- 
tum and classical theories. Extensive lists of the litera- 
ture on this subject can be found in reviews and books 



The ideas of Moyal were further developed in the late 
sixties in the works of Cahill and Glauber [jl0| and Agar- 
wal and Wolf pd| . As mentioned, the Wigner function 
is related to the Weyl (symmetric) ordering of the posi- 
tion and momentum operators q and p or, equivalently, 
of the bosonic annihilation and creation operators a and 

. However, there exist other possibilities of ordering. In 
particular, it was shown fin] that the Glauber-Sudarshan 



P function (ljJT^l is associated with the normal order- 
ing and the Husimi Q function |l4j| with the antinormal 
ordering of a and . Moreover, a whole family of s- 
parameterized functions can be defined on the complex 
plane which is equivalent to the q-p flat phase space. The 
index s is related to the corresponding ordering proce- 
dure of a and a'; the values +1, 0, and —1 of s corre- 
spond to the P, W, and Q functions, respectively. These 
phase-space functions are referred to as quasiprobabil- 
ity distributions (QPDs), as they play in the Moyal for- 
mulation of quantum mechanics a role similar to that 
of genuine probability distributions in classical statisti- 
cal mechanics. Various QPDs has been extensively used 
in many quantum-optical applications [ p)|Jil| ]. Most re- 
cently, there is great interest in the s-parameterized dis- 
tributions because of their role in modern schemes for 
measuring the quantum state of the radiation field [jlTj . 

The mathematical framework and the conceptual back- 
ground of the Moyal quantization have been essentially 
enlarged and generalized in two important papers by 
Bayen et al. |18|. Specifically, it was shown that noncom- 
mutative deformations of the algebra of classical phase- 
space functions (defined by the ordinary multiplication) 
give rise to operator algebras of quantum mechanics. 
This fact means that introducing noncommutative sym- 
bol calculus based on the so-called twisted product (also 
known as the star or Moyal product) , one obtains a com- 
pletely autonomous reformulation of quantum mechanics 
in terms of phase-space functions instead of Hilbert-space 
states and operators. This program of "quantization by 
deformation" has been developed in a number of works 

For a long time applications of the Moyal formulation 
were restricted to description of systems like a spinless 
non-relativistic quantum particle or a mode of the quan- 
tized radiation field (modeled by a quantum harmonic 
oscillator), i.e., to the case of the fiat phase space. Thcrc- 
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fore, an important problem is the generalization of the 
standard Moyal quantization for quantum systems pos- 
sessing an intrinsic group of symmetries, with the phase 
space being a homogeneous manifold on which the group 
of transformations acts transitively J2(],^4| . It has been 
recently understood that this problem can be solved us- 
ing the Stratonovich-Weyl (SW) correspondence. The 
idea of the SW correspondence is that the linear bijec- 
tive mapping between operators on the Hilbert space and 
functions on the phase space can be implemented by a 
kernel which satisfies a number of physically sensible pos- 
tulates, with covariance and traciality being the two most 
important ones. This idea first appeared in a paper by 
Stratonovich |25| in 1956, but it was almost forgotten 
for decades. The SW correspondence, that has been 
restated some years ago by Gracia-Bondia and Varilly 
|^6],^7|, has given a new impulse to the phase-space for- 
mulation of quantum theory. The SW method of the 
Moyal quantization has been applied to a number of im- 
portant situations: a nonrelativistic fr ee p article with 
spin, using the extended Galilei group [E6|; a relativis- 
tic free particle with spin, using the Poincare group p8[ ; 
the spin, using the SU(2) group |27]]; compact semisim- 
ple Lie groups p9| ; one- and two-dimensional kinemati- 
cal groups ]30| |33|] ; the two-dimensional Euclide an g roup 
[fjojpif ; systems of identical quantum particles [p5| . For 
a review of basic results see Ref. M . 

Notwithstanding the success of the SW method in the 
Moyal quantization of many important physical systems, 
the theory suffered from a serious problem. Specifically, 
it was the absence of a simple and effective method for the 
construction of the SW kernel which should implement 
the mapping between Hilbert-space operators and phase- 
space functions. The construction procedures for the SW 
kernels, considered during the last decade (see, e.g., Ref. 
||), did not guarantee that the kernel will satisfy all the 
SW postulates. Only very recently a general algorithm 
for constructing the SW kernel for quantum systems pos- 
sessing Lie-group symmetries was proposed | fj6| . It has 
been shown that the constructed kernel explicitly sat- 
isfies all the desired properties (the SW postulates) and 
that in the particular cases of the Heisenberg-Weyl group 
and SU(2) our general expression reduces to the known 
results. 

In the present paper we essentially extend the results 
of Ref. |36| and present a self-consistent theory of the 
SW method for the phase-space formulation of quan- 
tum mechanics. This theory makes use of the concept 
of generalized coherent states and of some basic ideas 
of harmonic analysis. Like the Cahill-Glauber formal- 
ism for the Heisenberg-Weyl group, we construct the s- 
parameterized family of functions on the phase space of 
a quantum system whose dynamical symmetry group is 
an arbitrary (finite-dimensional) Lie group. Accordingly, 
we introduce s-generalized versions of the traciality con- 
dition and the twisted product. The developed phase- 
space formulation is used for a general group-theoretic 
description of the quantum state reconstruction method. 



This description can be useful not only for measurements 
of quantum states but also in the field of signal process- 
ing. 



II. BASICS OF MOYAL QUANTIZATION 

A. Generalized coherent states and the definition of 
quantum phase space 

Given a specific physical system, the first thing one 
needs to do for the Moyal quantization (i.e., for con- 
structing phase-space functions) is to determine what is 
the related phase space. This can often be done by anal- 
ogy with the corresponding classical problem, thereby 
providing a direct route for the quantum-classical corre- 
spondence. From the technical point of view, the phase 
space can be conveniently determined using the concept 
of coherent states |37|]. The coherent-state approach is 
not just a convenient mathematical tool, but it also helps 
to understand how physical properties of the system are 
reflected by the geometrical structure of the related phase 
space. It is possible to say that the concept of coher- 
ent states constitutes a bridge between the Moyal phase- 
space quantization and the Berezin geometric quantiza- 
tion §| . 

Let G be a Lie group (connected and simply connected, 
with finite dimension n), which is the dynamical symme- 
try group of a given quantum system. Let T be a uni- 
tary irreducible representation of G acting on the Hilbert 
space H.. By choosing a fixed normalized reference state 
IV-o) £ "H, one can define the system of coherent states 



\ip g ) = r( ff )|Vo>, 



geG. 



(2.1) 



The isotropy subgroup H C G consists of all the group 
elements h that leave the reference state invariant up to 
a phase factor, 



T(h)\ipo 



|e^W| = l, heH. (2.2) 



For every element g € G, there is a decomposition of g 
into a product of two group elements, one in H and the 
other in the coset space X = G/H, 



Qh, 



geG, heH, neX. 



(2.3) 



It is clear that group elements g and g' with different 
h and h' but with the same Q produce coherent states 



which differ only by a phase factor: \ip g ) 



J5 



where 5 = 4>{h) — 4>{h'). Therefore a coherent state 
\Q) = \ipn) is determined by a point Q = Q(g) in the 
coset space X . A very important property is the identity 
resolution in terms of the coherent states: 



d/i(fi)|fi)(fi| 



(2.4) 



x 
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where dfi(fl) is the invariant integration measure on X, 
the integration is over the whole manifold X, and / is 
the identity operator on TL. The natural action of G on 
X will be denoted by g • Q,. 

An important class of coherent-state systems corre- 
sponds to the coset spaces X = G/H which are homo- 
geneous Kahlerian manifolds. Then X can be consid- 
ered as the phase space of a classical dynamical system, 
and the mapping £1 — > \Q) (p,\ is the geometric quantiza- 
tion for this system [|38| . The standard (or maximum- 
symmetry) systems of the coherent states correspond to 
the cases when an 'extreme' state of the representation 
Hilbcrt space (e.g., the vacuum state of an oscillator or 
the lowest /highest spin state) is chosen as the reference 
state. This choice of the reference state leads to systems 
consisting of states with properties "closest to those of 
classical states" J37]|3j| . In what follows we will consider 
the coherent states of maximal symmetry and assume 
that the phase space of the quantum system is a homoge- 
neous Kahlerian manifold X = G/H, each point of which 
corresponds to a coherent state |f2). In particular, the 
Glauber coherent states of the Heisenberg-Weyl group 
H3 are defined on the complex plane C = H3/U(l), and 
the spin coherent states are defined on the unit sphere 
§ 2 = SU(2)/U(1). In the more rigorous mathematical 
language of Kirillov's theory [^(J, the phase space X is 
defined as the coadjoint orbit associated with the uni- 
tary irreducible representation T of the group G on the 
Hilbert space TL. 



B. The Stratonovich-Weyl correspondence 

Once the phase space of a quantum system is deter- 
mined, the Moyal quantization proceeds in the following 
way. Let A be an operator on TL. Then A can be mapped 
by a family of functions F A (£l; s) onto the phase space 
X (the index s labels functions in the family). If A is 
the density matrix p of a quantum system, the corre- 
sponding phase-space functions F p (Q;s) = P(Cl;s) are 
called QPDs. Of course, the phase-space formulation 
of the quantum theory for a given physical system can 
be successful only if the functions F A (fl; s) possess some 
physically motivated properties. These properties were 
formulated by Stratonovich p5| and are referred to as 
the SW correspondence: 

(0) Linearity: A — > F A (Cl;s) is one-to-one linear map. 

(1) Reality: 



F Ai (n;s) = [F A (n;s)}*. 



(2.5a) 



(ii) Standardization: 

/ d^(fl)F A {fl; s) = Tr A. (2.5b) 
Jx 



(iii) Covariance: 

F A(g) (n-s) = F A {g-n- S ), (2.5c) 

where A{g)= Tig- 1 ) AT (g). 

(iv) Traciality: 

/ d(J,{Q)F A (Q;s)F B (Cl;-s) =Tr(AB). (2.5d) 
Jx 



If the function F A (Cl; s) satisfies the SW correspondence, 
it is called the SW symbol of the operator A. 

The above conditions have a clear physical meaning. 
The linearity and the traciality conditions are related to 
the statistical interpretation of the theory. If B is the 
density matrix (the stat e op erator) of a system, then 
the traciality condition ( 2.5oj ) assures that the statisti- 
cal average of the phase-space distribution F A coincides 
with the quantum expectation value of the operator A. 
O'Connell and Wigner have shown that the tracial- 
ity condition for density matrices of a spinless quantum 
particle (there it appears as an overlap relation) is nec- 
essary for the uniqueness of the definition of the Wigner 
function. It has been also shown J2i| that the traciality 
condition is necessary for the uniqueness of the defini- 
tion of the symbol calculus (twisted or "star" products) 
of the phase-space functions and for the validity of t he rc - 
lated non-commutative Fourier analysis. Equation (2.5d) 
is actua lly a g eneralization of the usual traciality condi- 
tion J2^,^,^9|, as it holds for any s and not only for the 
Wigner case s = 0. The reality condition (2.5a) means 
that i f A is self-adjoint, then F A (£l; s) is real. The condi- 
tion (2.5b) is a natural normalization, which means that 
the image of the identity operator I is t he constant func- 
tion 1. The covariance condition ( 2.5c ) means that the 
phase-space formulation must explicitly express the sym- 
metry of the system. 

The linearity is taken into account if we implement the 
map A — > F A (Sl; s) by the generalized Weyl rule 



F A (0;s)=Tr[AA(Q;s)] ) 



(2.6) 



where {A(r2;s)} is a family (labeled by s) of operator- 
valued functions on the phase space X. These operators 
are referred to as the SW kernels. The generalized tra- 



ciality condition (2.5d) is taken into account if we define 
the inverse of the generalized Weyl rule (|2.6|) as 



A= [ dfi(n)F A (n;s)A(n;-s). (2.7) 
Jx 



Now, the conditions (2.5a)-(2.5c) of the SW correspon- 
dence for F A (ri;s) can be translated into the following 
conditions on the SW kernel A(f2; s): 

(i) A(Q;s) = [A(f7;s)] t VSlel 



(ii) / d/i(0)A(fi;s) = /. 



A" 



(iii) A(g-n ]S ) = T(g)A(n; S )T(g- 1 ). 



(2.8a) 
(2.8b) 
(2.8c) 
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Substituting the inverted maps ( |2.7|) fo r A and B into 
the generalized traciality condition (|2.5d| ) , we obtain the 
relation between functions with different values of the 
index s: 

F A (Cl;s)= [ d(Ji(Sl')K SiS >(Sl,Sl')F A (Sl'; S '), (2.9) 



x 



K a>s , {SI, Si') = Tr [A(Sl; s)A(Q'; -s% 



(2.10) 



If we take s = s' in Eq. (|j|) and take into account the 
arbitrariness of A, we obtain the following relation 

A(fi;s)= / ^(f2')^(^') A (^'; s )> (2.11) 

where the function 

K(Sl, Si') = Tr [A(Sl; s)A(Sl'; -s)] (2.12) 
behaves like the delta function on the manifold X. 



III. CONSTRUCTION OF THE 
STRATONOVICH-WEYL KERNEL 

It is clear that the Moyal quantization for a physical 
system is accomplished by constructing the SW kernel 
A(Sl; s) that satisfies the SW postulates. Although the 
form of the SW kernel has been known for many systems, 
a general construction method was not known. A pro- 
cedure that was applied in many works p|,p0|-p2]| is as 
follows. An arbitrary point SIq 6 X is fixed and then an 
Ansatz is made for a self-adjoint operator A(ilo) (usu- 
ally only the case s = was c onsid ered) that satisfies 
the standardization condition ( 2.8b| ) and the following 
property: 



A(fio) = T( 7 )A(ft )T( 7 - 



V 7 e H Qo 



(3.1) 



where Hfi = {7 6 G | 7 • SIq = SIq} is the isotropy sub- 
group for Sl . For any SI E X there exists g E G such 
that g ■ SIq = SI, and then the SW kernel is defined by 



A(Q) = A(g • Sl ) = T{g)A{Sl )T{g- 1 ). 



(3.2) 



This kernel automatically satisfies the covariance condi- 



tion (2.8c), but the problem is that the traciality is not 
guaranteed. Of course, in the described procedure the 
form of the kernel depends on the Ansatz and often no 
kernel satisfying the traciality condition is found. 

We propose here a simple and general algorithm for 
constructing the SW kernels (the whole s-parameterized 
family) which explicitly satisfy all the SW postulates, 
including both the covariance and the traciality. Our 
method makes use of Perelomov's concept of coherent 
states and of only some basic ideas from harmonic anal- 
ysis. Hopefully, the simplicity and generality of our 
method will draw more attention to the ideas of the 
phase-space quantization. 



A. Necessary instruments: harmonic functions, 
invariant coefficients, and tensor operators 

Our problem is to find the explicit form of the SW ker- 
nel A (Sl; s) that satisfies the conditions (2.8a)-(2.8c) and 
( |2.1l| ). In order to accomplish this task, we need three 
basic ingredients: harmonic functions, invariant coeffi- 
cients, and tensor operators. The coherent states serve 
here as the glue that binds them together. 

We start by considering the Hilbert space L 2 (X, fi) of 
square-integrable functions u(Sl) on X with the invariant 
measure d\x. The representation T of the Lie group G on 
L 2 (X, /i) is defined as 



T(g)u(CL) = uig- 1 ■ Q). 



(3.3) 



The eigenf unctions Y V (SV) of the Laplace-Beltrami opera- 
tor [|42| form a complete orthonormal basis in L 2 (X, /i): 



5^y*(n)n(n / ) = *(fi-fi / ) 

dn(Sl)Y*(Sl)Y v ,(Sl) =<W 



(3.4a) 
(3.4b) 



A" 



The functions Y U (Q) are called the harmonic functions, 
and 5(n — SI') is the delta function in X with respect to 
the measure dfi. Note that the index v is multiple; it 
has one discrete part, while the other part is discrete for 
compact manifolds and continuous for noncompact man- 
ifolds. In the latter case the summation over v includes 
an integration with the Plancherel measure dp(v) and the 
symbol 8 vu i includes some Dirac delta functions (for more 
details see Ref. [^||). For conciseness, we omit these de- 
tails in our formulas. The harmonic functions Y v {Sl) are 
linear combinations of matrix elements T uu '{g). There- 
fore, the transformation rule for the harmonic functions 
is g§ 

T{g)Y v {Sl) = Y v (g- X ■ SI) =Y J T V ' V {g)Y v ,{St). (3.5) 



It should be understood that the summation in Eq. ( |3.5| ) 
is only on the part of v that labels functions within an 
irreducible subspace. 

Next, we once again use the coherent states, in order 
to introduce the concept of invariant coefficients. The 
positive- valued function |(f2|f2')| 2 is symmetric in SI and 
SI'. Therefore, its expansion in the orthonormal basis 
must be of the form 



\{SI\SI')\ 2 = J2t»Y u *(SI)Y v (SI>) 

= Y,t v y;{si')y v (si), 



(3.6) 



where r„ are real positive coefficients. Using the invari- 
ance (f2|f2'} = (g ■ Sl\g ■ SI') and the unitarity of the rep- 
resentation T, we obtain 
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|(Q|n')l 2 = E r - y ^-^(5-^) 

= J2 Y u*(^)J2 T ^'(9)Yu(9 ■ Of). (3.7) 



In order to satisfy this equality, the coefficients t v must 
be invariant under the index transformation of equation 
( |3 . 5[ ) : r„ = t v i . This means that t v do not depend on 
the part of v which labels functions within an irreducible 
subspace. Since the Laplace-Beltrami operator is self- 
adjoint, one finds that 



Y:(Q) = e i ^Y„(Q), 



(3.8) 



where Yp(O) is another harmonic function, with the same 
eigenvalue as Y V (Q). Since |(f2|f2'}| 2 is real, the coeffi- 
cients r„ must be inv ariant under the index transforma- 
tion of equation (|3.8| ): t v — t v . 

Next we use the coherent states, harmonic functions, 
and invariant coefficients for defining the set of operators 
{D v } on H: 



~ 1/2 / dn(tt)Y v (n)\n)(n\. 

J X 



(3.9) 



Using the e xpression (3.6) and the orthonormality re- 
lation ( |3.4b| ) for the harmonic functions, we obtain the 
orthonormality condition for the operators D u : 



Tr(D v Dl) 



(3.10) 



— 1/2 

Note tha t the factor r„ in front of the integral in 
Eq. (|3.9[) serves just for the proper normalization. Us- 
ing (|3.6j), we also obtain the relation 



T- l ' 2 {n\D v \n) = Y„(n). 



(3.11) 



The invariance of the coefficients t v implies that D v are 
the tensor operators whose transformation rule is the 
same as for the harmonic functions Y„(Q): 



T{g)D u T(g- 1 )=Y,T v > v {g)D p 



(3-12) 



A useful property of the tensor operators is that any op- 
erator A on TL can be expanded in the orthonormal basis 

A = Y J ^{ADt)D u . (3.13) 



B. Explicit form of the kernel 

Using the above preliminary results, we are able to find 
the SW kernel A(fi;s) with all the desired properties. 
Specifically, let us define 



A(n;s) = Y,f(s;r„)Y v *(n)D v . 



(3.14) 



We will s how that the construction of the generalized ker- 
nel (|3.14 ) satisfies the SW correspondence. In equation 
( ^.14| ) f(s; t v ) is a function of t v and of the index s. We 
assume that / possesses the invariance properties of r v . 

Using the inv ariance of r„ under the index transforma- 
tion of Eq. ( |3.8| ), we see that the reality condition (2.8a) 
is satisfied if f[s; t„) is a real-valued function. Therefore, 
we can consider only real values of the index s. 

Next we consider t he sta ndardization condition ( 2.8lj ) . 
Using the definition ( |3.14| ), we obtain 

f dM(fi)A(fi;s) =J2f(s; Tv )D v [ dn(Cl)X?(n), 
Jx v J x 

(3.15) 

while Eq. (|3.13D can be used to write 

v v J X 

(3.16) 

The s tanda rdizat ion c ondition is satisfied if the expres- 
sions ( 3.1 5|) and ( 3.1 6| ) are equal. Using the identity res- 
olution (2.4) and Eq. (3.6), we can write 



/\ |2 



i = (n\n) = / dn(n')\(n\n'}\ 

Jx 



= J2^ Y u(0) I dn(SV)Y v (tf). (3.17) 

Multiplying the left and right sides of this equation by 
Yy'(f2) and integrating over dfi(£l), we obtain 



dfj,(Cl)Y u (n) = t u I dfi{0)Y v (Cl). (3.18) 
x Jx 

Since r„ is not identically 1 , this relation can be satisfied 
only if there exists some i/q such that r„ = 1 and 



dfx(0)Y v (Q) oc 8 UV 



(3.19) 



A" 



(As was already mentioned, for noncompact manifolds 
the symbol 5 vv i actually includes some Dirac delta func- 
tions.) It can be easily seen from Eqs. ( 3.15| ), ( [3.16 ), and 
(|3.19|) that the standardization condition is satisfied if 

f(s',T V0 ) 



.-1/2 



/(s;l) = l, Vs. 



(3.20) 



The covariance condition (2.8c) can be rewritten as 
Y,f{s;r«)D v Y:{g-n) 

= J2f&Tu)T(g)D v T(g- 1 )Y;(n). (3.21) 
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Using the transformation rules (3.5) and (3.12), Eq. 
( |3.2l[) can be transformed into 



^^/( s; r„)T„v(#,y;(^). 



(3.22) 



Cha nging the summation indexes v «-► v' on either side of 
Eq. ( 3.22| ), we immediately see that the covariance con- 
dition is satisfied by the virtue of the inva riancc of r v 
under the index transformation of E qs. (|3.5| ) and ( |3.12| ). 

In order to satisfy the r elation ( 2.11| ), the function 
K(Cl,D,') of equation (2.12) must be the delta function 
in X with respect to the measure dfj,, 

K(n,n') = ^2Y*(n)Y v (n') =s(n-n'). (3.23) 



This result is valid if 

f(s;r v )f(-s;T u ) = 1. 



(3.24) 



This property is satisfied only by the exponential func- 
tion of s, i.e., 



/(*;<> = [/«)]* 



(3.25) 



Note that the standardization condition ( [3.2C ) then reads 
/(1) = 1- 

The exact form of the function f(r u ) can be determined 
if we define |43| for s = — 1 



A(fi ; -i) = \n)(n\. 

Then we obtain 

\n)(n\ = J2lf(T»T 1 Y v *mDu 



(3.26) 



(3.27) 



Multiplying both sides of this equation by Y v /(y£) and 

— 1/2 

integrating over e?/i(Sl), we find f(r v ) = t„ , i.e., 



f(s;T v ) 



s/2 



(3.28) 



Obviously, the standardization condition /(l) = 1 is sat- 
isfied. Finally, we obtain 



V 



(3.29) 



It is evident that this kernel is completely determined by 
the harmonic functions on the corresponding manifold 
and by the coherent states w hich f orm this manifold. We 
will see that the SW kernel ( 3.29] ) is a generalization of 
the Cahill-Glauber kernel for a harmonic oscillator [ ^0|pT| 
and of the Agarwal kernel for spin W4| . 



IV. PHASE-SPACE FUNCTIONS AND THE 
SYMBOL CALCULUS 

A. Types of phase-space function 

As the explicit form of the SW kernels is known, we 
can write the SW symbols on the phase space as 

v 

= £r-*/VU;(fi), (4.1) 

V 

where we have defined 

A v = Tr (ADD, X, = Tr (AD V ). (4.2) 

For a self-adjoint operator A, we get A v — A* v . It can 

be ea sily verified that substituting exp ressions (4.1) and 
(HH) into the inverse Weyl rule ( |2.7[ ), one indeed ob- 
tains A = ~^2 V A V D V . We a lso note that the function 
K s-S i (SI, Q') of equation (|2.10) is given by 



K s , s ,(n,n') = ^r-( s - s ')/ 2 y;(si)y;(si'), (4.3) 



and it clearly satisfies Eq. ( |2.9| ) which connects the func- 
tions with different values of the index s. In general, let 
F(fl) and if(Sl) be two phase-space functions such that 

2T(fi) = ^ j H v Y v {Sl). (4.5) 

V 

Then they are related through the transformation 

f(q,)= I dn(n')K FH {n,n')H(n'), (4.6) 

J X 

K FH (n,n>) = £ ^Y v {n)YZ{n'). (4.7) 



Let {| </>„)} be a complete orthonormal basis in the 
Hilbert space TC. Using the generalized Weyl rule fl2.6j ) 
for the operator A — \4> n ){4>m\, we obtain 



F A (Q;s) 



t |A(fi; s) \<f> n ) =A mn (n;s). (4.8) 



Using Eq. (|3.29|) , we find 

A TOn (0;s) = J2 T- 8 l 2 {<f> m \Dt \<f> n )Y v (n). (4.9) 



The s tandardization and traciality conditions (2.5b) and 
( |2.5d| ) can be used to show that 
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d^(Cl)A mn (Cl; s) = s n 



/ dfj,(Cl)A mn (Cl; s)A k [(tt; -s) = S m i5 n 
Jx 



(4.10) 
(4.11) 



w A (ci)*w B {ci) = w AB (ci) 



(4.18) 



The functions A mn (f2; s) form a useful orthonormal basis 
in L 2 {X, n). 

The SW symbols obtained for some special values of 
s are frequently used in numerous applications. In par- 
ticular, for s = — 1, we obtain the Q function (Berczin's 
covariant symbol p8[): 



Q A (si) = F A (ii;-i) = (cl\a\cl). 



(4.12) 



Equa tion ( 4.12 ) can b e easily obtained by recalling [see 
Eqs. Q3.26D and (|3.27|)] that 



A(n ; -i) = |n)<n| =Y,tI ,2 y:{V)d v . 



(4.13) 



For s = 1, we obtain the P function (Berezin's contravari- 
ant symbol ||): 

P A {Cl) = F A (Cl; 1) = ]T T^' 2 A v Y v (Ci), (4.14) 
whose defining property is 

a = [ dfi(n)p A {n)\n)(n\. (4.15) 

Jx 

The functions P and Q are counterparts in the tracial- 
ity condition (2.5d). Perhaps the most important SW 
symbol corresponds to s — 0, because this function is 
"self-conjugate" in the sense that i t is t he counterpart 
of itself in the traciality condition (2.5d). It is natural 
to call the function with s = the generalized Wigner 
function: 

W A {n) = F A (Q;0) = J2A v Y v (Q). (4.16) 

The corresponding SW kernel is 

A(ft;0) = A w (f>) = ]Ty;(f!)A,. (4.17) 



B. The generalized twisted product 

The phase-space formulation of quantum mechanics 
can be made completely autonomous if one introduces 
a symbol calculus for the functions on the phase space, 
which replaces the usual manipulations with operators 
on the Hilbert space. This symbol calculus is based on 
the so-called twisted product (or Moyal product) which 
corresponds to the usual product of operators |ll|^6|,^7j . 

Let us first consider the case of the Wigner function 
(s = 0). The twisted product of two functions is denoted 
by W A * Wb and is determined by the condition 



for an y two operators A and B. Note that the condition 
(4.18) assures the associativity of the twisted product. 
On the other hand, this product is, in general, noncom- 
mutative. In this way the algebra of operators is mapped 
onto the algebra of phase-space functions. If one starts 
from a classical phase-space description, the introduction 
of the twisted product can be viewed as the quantization 
realized by a deformation of the algebra of functions [18| . 



Using the Weyl rule (^6) and its inverse (2/7), we ob- 
tain 

W AB {n) = Tr[A w (Q)AB] 
= Tr 



A w (fl) I dn(n')W{Q.')A w (n') 

X 



x / d^{fi")w(n')A w (n") 

Jx 



(4.19) 



Introducing the function (trikerncl) 

L(Sl,Q',Q") = Tr[A w (Cl)A w (fl')A w {fl")}, (4.20) 
we obtain the following definition of the twisted product: 

(w A *w B ){n) = [ [ d^(ci')dfj,(Q")L(Q,n',n") 

Jx Jx 

xW A (n')W B (£l"). (4.21) 

The so-called Moyal bracket is defined as 

[W A , W b ]m = ~i{W A *W B -W B * W a ). (4.22) 

The twisted product can be easily generalized for ar- 
bitrary values of s. The s-parameterized twisted prod- 
uct (F A * Fb)(Q;s) of any two functions F A (Q,;s') and 
Fb(Q; s") is once again determined by the condition 



F A (n ; s ') * F B (n ; s ") = F AB (n- «). 



(4.23) 



Analogously to the Wigner function case, this leads to 
the definition 

(F A *F B )(n;s) = [ [ dti(ti')dn(n")L SjS , iS „(n,n',n") 

JxJx 

xF A (n'; S ')F B (n"; S "), (4.24) 
where the generalized trikernel is given by 

L S:S , tS „ (O, CI', CI") = Tr [A(fi; s)A(ft'; -s')A(0"; -s")] 
= A mn (Cl; S )A nk (Cl';~s')A km (Cl"- 1 -s"). (4.25) 



Using the standardization condition (2.8b) and the def- 
inition (2.1C), we obtain 
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/ d/i(fi)L s , s /, a »(0,fi',fi") = Tr[A(Q';-s')A(Q";-s")] 
Jx 

= K_ s ,. s „{Q',n"). (4.26) 

This result together with the relation can be used 
to obtain the so-called tracial identity for the generalized 
twisted product, 



dii(Sl)(F A *F B )(n;8) 



x 



[ d^n)F A (Q;s')F B (n;-s'), (4.27) 
Jx 



which holds for any s and s'. Equation ( ggg ) is the i ze d in the following way 



field (modeled by a quantum harmonic oscillator). The 
Wigner function and the Moyal quantization g] were 
originally introduced for such systems. The kernel im- 
plementing the mapping between Hilbcrt-space opera- 
tors and s-parameterized families of phase-space func- 
tions (the SW kernel in our notation) for H3 was intro- 
duced by Cahill and Glauber |l0| . The generalization of 
the formalism to the many-dimensional case is straight- 
forward (see, e.g., Ref. ||). 

The nilpotent Lie algebra of H3 is spanned by the basis 
{a, a^, /}, where a and are the boson annihilation and 
creation operators, satisfying the canonical commutation 
relation, [a, a^] = /. Group elements can be parameter- 



phase-space version of the tracial identity for the opera- 
tors, 



Tr (AB) = A V B V . 



(4.28) 



Usi ng the covariance condition ( p.8q) and the definition 
(4.25), we find the invariance property of the trikernel 



L S:S ,, S „ (g-n,g-n',g- O") = L S:S ^ S „ (0, O', O"). (4.29) 

This property implies the equivariance of the twisted 
product: 

(F A * F B y(n; s) = s') * Fg(fi; s"), (4.30) 

where 

F g A (Q;s) = F A (g- 1 -n;s). (4.31) 



V. EXAMPLES 

The general formalism presented above can be under- 
stood much better by illustrating it with a number of 
simple examples. We will consider two simple physi- 
cal systems: a (nonrelativistic) spinless quantum parti- 
cle and spin, whose dynamical symmetry groups are the 
Heisenberg-Weyl group H3 and SU(2), respectively. It 
should be emphasized that the SW kernels for these ba- 
sic systems have been known for a long time ^5| , so the 
novelty here is not the result itself but the method of 
derivation. Our aim is to demonstrate how the general 
algorithm works by applying it to a number of relatively 
simple and well-known problems. We will show that by 
identifying harmonic functions, invariant coefficients, and 
tensor operators for a given system, one can readily de- 
rive the explicit form of the SW kernel. 



A. The Heisenberg-Weyl group 

First, we consider the Heisenberg-Weyl group H 3 which 
is the dynamical symmetry group for a spinless quan- 
tum particle and for a mode of the quantized radiation 



where 7 £ C and (p £ R. 

The phase space is the complex plane C : 
and the (Glauber) coherent states are 



|fi> = |a> =D(a)|0), aeC, 



where 



(5.1) 
Ha/U(l), 
(5.2) 

(5.3) 



D(ot) — exp(aa^ — a*a) 

is the displacement operator. The invariant measure is 

dn(n) = Tr~ 1 d 2 a, (5.4) 

and the corresponding delta function is 

<5(n-n') =<irSW {a -a'). (5.5) 

The harmonic functions on C are the exponentials: 

Y v (tt) = Y 6 (a) = T(e, a) = exp(£a* - (5.6) 

Here v = £ 6 C with the Plancherel measure given by 
dp{v) = -K- X d 2 i and with 6 v y = nS^ - £). Note that 
for the Heisenberg-Weyl group both the phase-space co- 
ordinate = a and the index v = £ are complex num- 
bers, and the Plancherel measure is similar to the invari- 
ant measure on C. 

The invariant coefficients t v = r(£) can b e fou nd in 
the following way. In the present context Eq. (3.6) takes 
the form 



/\ |2 



\(a\a')\ 



£*(a— a )— — a )* 



Taking into account that the Fourier transform of a Gaus- 
sian function is once again a Gaussian, it is not difficult 
to obtain 

t(£)=cx P (-|£| 2 ). (5.7) 
Then we deduce that the tensor operator 
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D„ = £>(0 



d 2 a 



a)(a\ (5. 



is just the displacement operator £>(£) = e^ a ~^ a . The 
natural orthonormal basis in the Hilbert space is the Fock 
basis {|n)}, <va|n) = n\n) (n = 0, 1, 2, . . .). The matrix 
elements of the tensor operator are given by p7| 

(m\D(0\n) 



( y/rHJrrile-^ /2 £ m ~ n L%- n (\£\ 2 ), m>n 

\ y^!7de-l«l 2 /2(_^)»-™i«-»'(|^|2) i m < n , 

where L^(x) are the associ ated Laguerre polynomials. 
Using the parameterization (5.1) of group elements, one 
can easily find the transformation rule 



T(g)D(0T(g- 1 ) = D(7)D(^)D(— 7) 

= exp( 7 r- 7*0^(0- 



(5.9) 



Therefore, the index £ does not change under the group 
transformation, as D(£) and Y(£,a) are just multiplied 
by a phase factor. Correspondingly, there is no index 
transformation, induced by the action of group elements, 
to which r(£) should be invariant. On the other hand, 
a) — Y(— £, a), and r(£) is obviously invariant un- 
der the index transformation £ «-> — £. 

Finally, the harmonic functions Y(£, a), the invariant 
coefficients r(£), and the tensor oper ators 13(0 can be 
substituted into the general formula ( 3. 29] ). Then one 
obtains the SW kernel for the Heisenberg-Weyl group: 



A(a;s)= I ^X s l«l 2 /V* Q -« Q *e« at -«* Q , 

' 7T 



(5.10) 



which is exactly the kernel introduced by Cahill and 
Glauber @. 



B. The SU(2) group 

As another example, we consider SU(2) which is the 
dynamical symmetry group for the angular momentum 
or spin and for many other systems, for example, a 
collection of two-level atoms, the Stokes operators de- 
scribing the polarization of the quantized light field, 
two light modes with a fixed total photon number, etc. 
A number of authors have used different approaches 
to the construction of the Wigner function for spin 
@f||23H@IIIH. The explicit expressions for the 
Q, W, and P functions for arbitrary spin were first ob- 
tained by Agarwal f44||, w ho used the spin coherent-state 
representation [3^j53| , |54[ | and the Fano multipole opera- 
tors J55| . Varilly and Gracia-Bondfa (27) have shown that 
the spin coherent-state approach is equivalent to the for- 
malism based on the SW correspondence. 

The simple Lie algebra of SU(2) is spanned by the basis 

{J xi Jyi Jz\l 



(5.11) 



The unitary irreducible representations are labeled by 
the index j (j — 0,1/2,1,...), and the Hilbert space 
TLj is spanned by the orthonormal basis \j, ji) (ji — 
j,j — 1, . . . , —j). Group elements can be parameterized 
using the Euler angles a, /?, 7: 



(5.12) 



The phase space is the unit sphere S 2 = SU(2)/U(1), 
and each coherent state is characterized by the unit vec- 
tor 



n = (sin 6 cos 4>, sin 9 sin cf>, cos ( 



(5.13) 



Specifically, the coherent states |Q) = \j; n) are given by 
the action of the group element 

g{Q) = g(Q, 0) = e-^J'e-^Jy (5.14) 

on the highest-weight state 

\j;n) = \j;6,<j>}=g(d,ct>)\j,j) 



^ ( 2j ^ 2 



E 



3 + (j> 



(0/2) sin J '-"(0/2) 



xe-^|i,/x). 
The invariant measure is 



(5.15) 



dfjbtn) = 2j ' + 1 dn = sin dd dd>, (5.16) 

47T 47T 

and the corresponding delta function is 
4tt 



s(n - n') 



2j + l 

47T 

2j + l 



<5(n - n') 
5 (cos - cos 6' ) 5 (0 -ft). (5.17) 



The harmonic functions on S 2 are the familiar spherical 
harmonics: 



47T 
27Tl 



(5.18) 



In this context v is the double discrete index {I, m} with 
I = 0, 1, 2, . . . and m = I, I — 1, . . . , — /. The transforma- 
tion rule for the spherical harmonics reads 



(5.19) 



where 



ft 7) = {l,m'\g(a,(3 7l )\l,m) (5.20) 
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is the matrix representati on o f SU(2) elements and 
g(a,P,"f) is given by Eq. ( 5.12 ). Another property of 
the spherical harmonics is 



Y, 



irVTYi,- 



(5.21) 



The invariant coefficients can be found using the fol- 
lowing expansion |27[]: 



\{jMj^)\ = 



2 , 1 + n-n' 



2j 



2,7 



Q/4-1 

E JrrrrO*>i; «, j> 2 fl(n • n'), (5.22) 



2=0 



2j + l 



where Pi(x) are the Legendre polynomials and 

(ji,mi;j2,TO 2 |j,m) = C^ 2TO (5.23) 

are the Clebsch-Gordan coefficients. Using the addition 
formula for the spherical harmonics, 



2Z + 1 

4-7T 



P/(n-n')= E ^mW^mCn') 



(5.24) 



---i 



Eq. (5.22) can be rewritten as 



I0>|i,n')| 



4?r 



2j J 



J i=o m=-/ 

xy4(n)y Zm (n') 



(5.25) 



Comparing this result with the general formula ( |3.6| ), we 
readily find that the invariant coefficients are given by 



t v = n = (j, 



(2j + l)[(2j)!] 2 
(2j + i + l)!(2j-i)!' 



(5.26) 



Note that 77 = for I > 2j. The invariance of 77 is 
ensured by the fact that they are independent of m. 

The tensor operators for spin are the well-known Fano 
multipole operators pqj , which can be written in the form 



Di m — 



21 + 1 



2j 



- E (3,k;l,m\j,q}\j,q){j,k\. (5.27) 



k,q=-j 



Substituting expressi ons lyj), ( ]5~26| ), and ( |5~27| ) into 
the general formula ( 3.29 ), we find that the SW kernel 
for spin is given by 



4tt 



2j 



-E ./-./ :/ - "./•./ s 

2=0 

E Am^M), (5.28) 



which coincides for s = 0, ±1 with the results by Agarwal 
ffl and by Varilly and Gracia-Bondi'a p7|. 



VI. RECONSTRUCTION OF QUANTUM STATES 

A. Basic systems and methods 

A great amount of work has been devoted in the last 
few years to the problem of determining the quantum 
state from information obtained by a set of measure- 
ments performed on an ensemble of identically prepared 
systems. The task is to reconstruct the density matrix 
p which, according to the principles of quantum physics, 
contains all available information about the state of a 
system. Of course, the question arises which set of mea- 
surements provides information sufficient for the state re- 
construction. This question first appeared in early works 
by Fano ]5t| a nd Pauli |57j| and was discussed in a number 
of papers i|pq-B3j . 

Recently, significant theoretical and experimental 
progress has been achieved in the reconstruction of quan- 
tum states of the light field (see, e.g., a recent book p7|). 
One of the most successful reconstruction methods in this 
context is the optical homodyne tomography. A tomo- 
graphic approach to the Wigner function was discussed 
by Bertrand and Bertrand p4| and a quantum-optical 
scheme was proposed by Vogel and Risken |pq| . The 
reconstruction of quantum states of the light field by 
means of homodyne tomog raphy was realized in a series 
of intriguing experiments |p6| , p7| . Various methods for 
data analysis in optical homodyne tomography measure- 
ments were recently discussed [^-[72] . The tomographic 
schemes were also generalized for the reconstruction of 
the joint density matrix for two- and multi-mode opti- 
cal fields 1 73 - 77 1 . Among other approaches to the re- 



construction of quantum states of light we would like to 
mention the symplectic tomography J78| and the pho- 
ton counting method s ]79| ~ |8l|l (also known as the photon 
number tomography |82[). 

In the case of a single-mode microwave field inside a 
high-Q cavity, a direct measurement on the system it- 
self is impossible. Instead, one can probe the state of the 
intra-cavity field via the detection of atoms after their in- 
teraction with the field mode [|S3|-[S5f . Similar ideas were 
also applied to the reconstruction of the quantum mo- 
tional state of a laser-cooled ion trapped in a harmonic 
potential [R3J86H90], including a beautiful experimental 
realization |pl|. 

State reconstruction procedures were proposed for var- 
ious quantum systems, for example, one-dimensional 
wave packets ||9^,P3|1, harmonic and anharmonic molec- 
ular vibrations [p4|]95| , motional states of atom beams 
p6| , Bose-Einstcin condensates p7t , cyclotron states of 
a trapped electron [|98|| , atomic Rydberg wave func- 
tions p9fl , etc. State reconstruction methods for sys- 
tems with a finite-dimensional state space (e.g., for spin) 
were also discussed pl| , p2| , p9[p0 10C. 101]. Experimental 
reconstructions were also reported for electronic angular- 
momentum states of hydrogen 102] an d for vibrational 
quantum states of a diatomic molecule [ 103 1 . 
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B. Displaced projectors 



r n (a) =D(a)\n){n\D\a) 



(6.6) 



It turns out that the majority of schemes used for 
the reconstruction of quantum states are related to the 
phase-space formalism. Frequently, the Q function, the 
Wigner function, or other phase-space QPDs represent- 
ing the density matrix p of the system can be cither 
measured directly or deduced in some way from mea- 
sured data. In particular, in many proposed and realized 
schemes the measured quantity is the expectation value 



P«(A) = (T„(A)> = Tr [ P T U (X)} 
of a self-adjoint operator 

r u (A) = U{X)\u){u\U\\), 



(6.1) 



(6.2) 



which is a transformed projector on a quantum state \u). 
The unitary operator U (A) represents the corresponding 
transformation, and the measurements are made for a 
range of values of the transformation parameter A. 

We will distinguish here between two possibilities. If 
U{X) = T(f2) is the phase-space displacement operator 
which represents an element of X = G/H, with G being 
the dynamical symmetry group of given quantum system, 
we will call the observable r u (A) = r„(Q) the properly 
transformed projector (or the displaced projector). Oth- 
erwise r„(A) will be called the improperly transformed 
projector. 

In order to illustrate these definitions, let us consider a 
quantum harmonic oscillator which is the model system 
for a single mode of the quantized radiation field, a laser- 
cooled ion moving in a harmonic trap, or a harmonic vi- 
brational mode of a diatomic molecule. The correspond- 
ing symmetry group is the Heisenberg-Weyl group H3, 
and the phas e spa ce is the complex plane C = H3/U(l) 
(see section VA). In this context U(X) = D(a) is 
the Glauber displacement operator, and the expectation 
value of the displaced projector, 

Pu (a) = Tr [pT u (a)] = Tr [pD(a)p u D^ (a)], (6.3) 

is calle d the ope rational phase-space probability distri- 
bution 104 -10£]. Here, p is the density matrix of the 
quantum state of the system and p u is the density ma- 
trix (given by the projector \u)(u\ for a pure state) of the 
so-called "quantum ruler" state which characterizes the 
measurement device. For example, displacing the state 
of the oscillator, 



p^p{a) = D^(a)pD(a), 



a G 



(6.4) 



and measuring the probability of finding it in the Fock 
state |n), one obtains the operational phase-space prob- 
ability distribution, 



Pn(a) = (n\p(a)\n) = Tr [pT n (a)}. 
The displaced projector 



(6.5) 



is obtained for \u) = \n) being the Fock state. In par- 
ticular, measuring the probability to find the displaced 
oscillator in the ground state |0), one obtains the Husimi 
function Q(a) — (a\p\a). On the other hand, if one 
knows the functions p n (a) for all values of n, then the 
Wigner function can be built as Jl0| 



W(a)= 2 £(- !)>„(«) 



(6.7) 



71 = 



This formu la ca n be generalized for QPDs with other 
values of s [ 107 1 : 



9 1 — ' / 4_ 1 \ 

F p (cr,8)=P(ar,8) = -—Y,(j=i) ( 6 ' 8 ) 

71=0 ^ ' 

These methods for determining the Husimi function 
and the Wigner function (and thus reconstructing the 
quantum state of the system) were first discussed by 
Royer [^2) in 1985. Recently, such a scheme for measur- 
ing the Q function was proposed in the context of trapped 
ions [^Tj. Another method for the reconstruction of the 
motional state of a trapped ion, proposed and experi- 
mentally realized by the NIST group |)lj], employs the 
interaction between the vibrational mode of the ion and 
its internal electronic levels. The initial mot ional state 
is displaced in the phase space, as in Eq. d6.4[ ), and then 
the interaction with the two-level internal subsystem is 
induced for a time t. The population P^ (t, a) of the lower 
internal state | |) is measured for different values of dis- 
placement amplitude a and time t (this measurement can 
be made with great accuracy by monitoring the fluores- 
cence produced in driving a resonant dipole transition). 
If I J.) is the internal state at t = 0, the signal averaged 
over many measurements is 



n(*»«) = 5 



00 

1 + 2jf>„(a) cos(2il rhn+ it)e 

n=0 



-7nt 



, (6.9) 



where f2 n! „ + i are the Rabi frequencies and j n are the 
experimentally determined decay constants. This rela- 
tion allows to determine the populations p„(a) of the 
displaced motional eigenstates. As one can see from 
Eq. ( |6.8D , the functions p n {at) in their turn can be used to 
calculate the QPDs P(a;s) (e.g., the Wigner function). 
Alternatively, the density matrix in the Fock representa- 
tion can be deduced directly from p n (a) . 

In the optical domain, the function p n {a) can be de- 
termined in principle as the probability of recording n 
counts with an ideal photodetector exposed to the dis- 
placed light field. In practice, one could use the unbal- 
anced homodyning detection ]79|-|82]| , in which the signal 
field is mixed in a beam splitter with the local oscillator 
of coherent amplitude j3 and the photon statistics of the 
superimposed field is then counted by a photodetector 
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of quantum efficiency £. The resulting counting statis- 
tics is denoted by p n (a,rf), where a — —Rj3/T is the 
effective displacement amplitude (T and R are the trans- 
mission and reflection coefficients of the beam splitter) 
and 77 = C| T 2 is the over all q uantum efficiency. In this 
realistic situation formula ( |6.8[ ) should be replaced by the 
following result M: 



P(a;s) 



— Y 

n=0 



2 + 77(3-1) 

V(s - 1) 



p n (a,r)). (6.10) 



This method of state reconstruction is sometimes called 
the photon number tomography. 

As an example of measurements with improperly trans- 
formed projectors, we mention the optical homodyne to- 
mography 1 65 6^] in which one measures the probability 
distribution P(x, 9) for the rotated field quadrature 



x g = x cos 9 + p sin 9 = U(9)xU i (9) . 



(6.11) 



The field quadratures x and p can be viewed as the scaled 
position and momentum operators of the harmonic oscil- 
lator, with a = 2~ x l 2 {x + ip), and U(6) — exp(i9a^a) 
is the rotation operator (known in optics as the phase 
shifter) on the phase plane. U(9) represents an element 
of the SO(2)~U(l) subgroup of the oscillator group H4 
whose algebra is spanned by {/, a, , a^a}. The improp- 
erly transformed projector is given by 



T x {6) = U{9)\x)(x\U\9), 



(6.12) 



where \x) are the position eigenstates. The measured dis- 
tribution P(x, 9) can be used for determinin g th e Wigner 
function via the inverse Radon transform |34|-^6|]. Al- 
ternatively, the density matrix in some basis (e.g., in 
the Fock basis) can be deduced directly from P(x, 9) by 
averaging a set of pattern functions fl68H7l|. Another 
example of measurements with improperly transformed 
projectors is the symplectic tomography ifrsj , in which 
the phase-space rotation is accompanied by the squeez- 
ing transformation. 

In the case of measurements with improperly trans- 
formed projectors, a general group-theoretic approach is 
problematic, because the number of possible transforma- 
tions is very large and one should consider each situa- 
tion separately. On the other hand, the method of prop- 
erly transformed projectors works uniformly for physical 
systems with different symmetry groups. For example, 
in the case of the SU(2) symmet ry (e .g., spin, two-level 



atoms, etc.), proposals appeared [100,101] for measuring 
the Q function, 

Q(n) = {j,n\p\j,n) = Tr(p|j,n)(j,n|), (6.13) 

or, more generally, for measuring the probability 

p M (n)=Tr[ P r M (n)], (6.14) 
r M (n)= 5 (n)|i, M )(j )/J | 5 - 1 (n) (6.15) 



of finding the displaced system 

p(n) =g- 1 (n)pg(n), ne§ 2 



(6.16) 



in the state \j, p) . These ideas for spin are conceptu- 
ally very similar to the proposals in the context of op- 
tical fields or trapped ions. Therefore, it seems natural 
to apply the phase-space formalism developed above to 
the general group-theoretic description of the state recon- 
struction method based on the measurement of displaced 
projectors. 



C. General reconstruction formalism 

From the practical point of view, the reconstruction 
procedure consists of two steps. First, the system de- 
scribed by the density matrix p is displaced in the phase 
space: 

p->p(ty =T-\n)pT(Q), Slgl (6.17) 

The second step is the measurement of the probability to 
find the (displaced) system in a quantum state \u), 



Pu (n) = (u\ P (n)\u) 



(6.18) 



Repeating this procedure for a large number of phase- 
space points f2, one can in principle determine the func- 
tion p u (Q). 



1. More about displaced projectors 

The information contained in the function p u {Q) is 
enough for the reconstruction of the density matrix. It is 
convenient to analyze this problem with the help of the 
displaced projector, 



r„(n) = T{n)\u){u\T-\n), 



(6.19) 



whose expectatio n v alue gives the measured function 
Pu(f2), as in Eq. (|6.l|). The displaced projector satisfies 
a number of useful properties: 



(i) It is a self-adjoint operator, 

r£(n) = r„(n) vnex 



(6.20) 



Since p u (£l) is not only real but also non- negative 
(this is a probability), r u (0) is also a non- 
negatively defined operator. 

(ii) Provided that the state \u) is normalized, r u (f2) is 
a trace-class operator of trace one, and the follow- 
ing standardization condition holds, 

/ dp,(fl)T u (Q) = I. (6.21) 
Jx 



This implies the normalization of p u (Q), 



x 



dp{Q)p u (Q) = 1. 



(6.22) 
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(iii) The displaced projector is manifestly covariant, 

T(g)T u (Q)T(g-i) = T u (g - SI). (6.23) 

Consequently, if p u (Sl) corresponds to the initial 
density matrix p, the function p u {g • SI) will cor- 
respond to the transformed density matrix p(g) = 
T(g-i)pT(g). 

Denoting the density matrix of the quantum ruler state 
by p u (which is for a pure state), the operational 

phase-space probability distribution reads 



p u {Sl) = Tr[pT(Cl)p u T- 1 (n)}. 



(6.24) 



Using the inverse Weyl rule ( |2.7| ) for the density matrix 
p, we obtain 

p u (Sl) = / dp J (Sl')P(Sl';s)Ti[T- 1 (Sl)A(Sl';-s)T(Sl)p u }, 
Jx 

where P(Sl; s) = F p (Sl; s ) is t he SW symbol of p. Now, 
the covariance property ( 2.8c ) can be used to obtain the 
following expression: 



p u (Sl) = [ dn(SY)P(Q ■ Si'; s)P u (Sl'; 
Jx 



-«), (6-25) 



where P u (Sl;s) is the SW symbol of p u . Therefore, the 
operational phase-space probability distribution p u (Sl) is 
given by a convolution of the two QPDs representing the 
quantum state of the system and the quantum ruler state 
of the measurement apparatus. In the particular case of 
the Heisenb erg-W eyl group and for s = 0, t he ge neral 
expression (3.25) reduces to the known result ]105[ 



Pu(a) 



-W(a + a')W u {a'). 



(6.26) 



If the quantum ruler state \u) — \ipo) is the reference 
state of the coherent-state basis, then 



r^(n) = |n>(n| = A(i2;-i) 

is the SW kernel with s = — 1, and 

p^ (Q) = (Sl\p\Q) = Q p (Q) 



(6.27) 



(6.28) 



is the Q function. However, except for this coincidence, 
the displaced projectors are not the SW kernels, as they 
do not satisfy the traciality condition. On the other hand, 
the functions p u (Sl) differ from the majority of QPDs, as 
they are positive on the whole phase space (which re- 
flects the fact that they are measurable probabilities). 
Usually the state \u) is chosen to belong to some com- 
plete orthonormal basis {| (fin)} which consists of energy 
eigenstates of a natural Hamiltonian of the physical sys- 
tem (e.g., the Fock basis for a harmonic oscillator or J z 
eigenstates for spin). Then there exists the relation 

5>*»(fi) = l, (6-29) 

n 

which follows from the completeness of the basis. 



2. Entropy 



A useful quantity for analyzing statistical properties 
of the system (in particular, the amount of noise) is the 
entropy. A phase-space version of the entropy can be 
introduced in the following way, 



S u = -I dfj,(Sl)p u (Sl) tapu(fi). 
fx 



For \u) = \tfi ), Eq. (6.30) gives 



x 



dn{Sl)Q p (Sl) InQ p (n), 



(6.30) 



(6.31) 



which is a generalization of the Wehrl entropy [108] that 
was defined originally on the flat p hase space of the Wcyl- 
Heisenberg group. The entropy fl6.30p can be useful in 
the reconstruction procedure, as it is a sensitive measure 
of the noise added to the system during the displacement 
and detection processes. A si milar situation exists also in 
the field of signal processing [ 109 1 . There \u) represents 
the test signal and p u (Sl) is a distribution on the time- 
frequency space. One can produce various test signals \u) 
and compute the corresponding entropies S u . Choosing 
\u) that minimizes the entropy, one obtains the optimal 
form of pattern analysis (in particular, this method al- 
lows to achieve data compression). 



3. Harmonic expansions 



A useful expression for p u (Sl) can be derived in the 
following way. Using the expansion 



P 



(6.32) 



we obtain 

Ptt(fi) = ^'R u {u\T- 1 {Sl)D u T{Sl)\u) 

V 

= E^E T ^( n )HA>'M. (6.33) 



v v' 



Expanding p u (Sl) in the basis of harmonic functions, 

p u (Sl)=J2rt ) K v Y v (Sl) ) (6.34) 



(u) 

we identify the coefficients n\, by means of the following 
relation 



(6.35) 
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Formally, we can write 
~ T ^ 1/2 J 

xY v (n ■ n')\(u\n')\ 



/,:; r, - i i d^(si)dii{VL')Y*(si) 

X J X 

,/\|2 



(6.36) 



but actually Eq. ( 6.35 ) is more convenient for calculating 
the coefficients 



Equation (6.34) for the functions p„(fi) corresponds to 
the expansion 

r u (fi) = i^Y:{VL)D v = J2 K^Y v {VL)Dt (6.37) 

for the displaced projectors. It follows from the proper- 
ties of r u (fi) that the coefficients n„ are positive and 
possess the same invariance properties as r v . Using the 
general result ( |4.6| ), we obtain the relation between the 
functions p u (Q) and p v (£i), corresponding to different 
quantum ruler states \u) and \v), 

Pu {n)= [ dn(n')K uv (n,n')p v (n'), (6.38) 

Jx 

(«) 

K UV (W) =J2^Y u (a)Y*(Q'). (6.39) 



4- Deducing the density matrix and quasiprobabilities 

Knowledge of the phase-space function p u (Cl) allows to 
reconstruct the density matrix in a simple way: 



dfi(n) Pu (n)Y;(n). 



(6.40) 



x 



Formally, we can also represent the density matrix by 
means of an integral transform of the displaced projec- 
tor: 



P = I d/x(fi)r„(n)r„(fi). 

Ix 



(6.41) 



This relation gives the density matrix in terms of a phase- 
space function r u (Tl), and in this sense it is the inverse of 
Eq. (S.l). The function r u (fl) is defined by its harmonic 
expansion, 



t (Q) =J2[* { v ) r 1 K»Y u (n). 



(6.42) 



We also obtain the following relation between the func- 
tions 7\i(0) and p u (fl), 



p u (fl)= / dfi(n')r u (fl')Tr[r u (Q)r u (Q% (6.43) 
Jx 



Tr[r u (0)1^(0')] - |( lt |T- 1 (0)T(n')|^}| 2 

= £4 u >4«>y,,(n)y;(n') 



(6.44) 



Certainly, the most convenient way for deducing the den- 
sity matrix from the measured functi ons p lt (f2) is by cal- 
culating the coefficients 1Z U via Eq. ( |6T0l) . 

The measured functions p u {Q) can be used also for the 
reconstruction of various QPDs which represent the den- 
sity matrix in the phase- spa ce formulation. According to 
the general expression (fhl|), the QPDs for the density 
matrix p are given by the harmonic expansion 



F p (ft; s) = P(ft; s) = J2 t^^U^SI). (6.45) 



Therefor e, on e can just use the coefficients 1Z U calculated 
via Eq. ( 3.40 ). On the other hand, Eq. ([0]) can be used 
to obtain the relation between the QPDs P(f2; s) and the 
measured functions p u (fl): 



p(n-,s) = / dix(n')K- s (n,n')p a (n'), (6.46) 



Pu (n)= dfi(n')K+(n,n')P(n'; S ), (6.47) 

Jx 



where the transformation kernels are 



KtM n') = E ±X YuMYW). (6.48) 



It can be easily shown that K+ S (fl, Q') — P^ft -1 -fi'; — s) 
where P u (Sl; s) is the SW symb ol of \u){u\, so Eq. ( |6.47j ) 
is consistent with the relation ( |6.25| ). 

As was already mentioned, if the state \u) is the 
reference state |V>o) of the coherent-state basis, then 
T^ (n) — \Q,)(Cl\ and the measured function p^ (f2) coin- 
cides with the function Q p (ft) = P(ft;—l). Comparing 
the harmonic expansions ( |3~34| ) and ( ]6T45| ) for the case 
\u) = \ipo), we hud the following relation: 



Ma) _ T l/2 



(6.49) 



In t his ca se we also obtain that the function r^ (CI) of 
Eq. ( |6.42 ) is just the P function, 



r i , (Cl) = P p (Cl) = P(Cl;l). 



(6.50) 



where 



Note that in the case of the Heisenberg-Weyl group 
can also calculate the QPDs using the formula 
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5. Examples 

We see that the mathematical procedure of the recon- 
struction of the density matrix p and its QPDs P(Sl; s) 
from the measured probabilit y p u (Q) actually consists 
of the simple transformation ( |6.40 ). The mathematical 
tools one needs for this procedure are the harmonic func- 
tions Y„(Q) and the invariant coefficients t v and Kv . 
In what follows we compute the explicit form of Ku for 
simple but instructive examples of the Heisenberg-Weyl 
group (with \u) being a Fock state) and the SU(2) group 
(with \u) being a J z eigenstate). 

In the case of the Heisenberg-Weyl group, we consider 
the probability p n {a) to find the displaced initi al sta te in 
the Fock state \n) (n = 0,1,2,.. .). Then Eq. fl6.35|) can 
be rewritten in the form 

K (»)(£)y(£, a ) = (n\D(-a)D(£)D(a)\n). (6.51) 



Using Eq. 



, we obtain 



D(-a)D(£)D(a) =Y(£,a)D{t). 
Therefore, the k coefficients are given by 

«(">($ = (n\D(0\n) = exp(-±|£| 2 )L„(|£| 2 ). 



(6.52) 



(6.53) 



Of course, for n = one gets /c^(£) = [t(^)] 1 ^ 2 . 

In the case of the SU(2) group, we consider the proba- 
bility p^(9, 4>) to find the displaced initial state in t he J z 
eigenstate \j,n) (/X = j, j - 1, • • • , -j)- Then Eq. ( |6.35| ) 
takes the following form 



2j + l 

47T 



0',M|3 1 (O,0)Di m g(9,<p)\j,fi). 

(6.54) 



Using the paramete rizati on (5.14) for g(6,<p) and the 
transformation rule ( |5.19 ), we can write 



g- 1 (6,<P)D lm g(9,<p)= ]T X>$ m (0,M)Am'- (6-55) 



Since the matrix element of the tensor operator, 



(j,f4 D im'\j,v) = J^j^(j,fJ--,l,m'\j,fi), (6.56) 



vanishes unless ml — 0, Eq. ( 6.54 ) reads 

-z>$,(o,M)0\ttJ,o|j,/i>. 



21 + !^(0 



Taking into account the fact that 



V i i ) m (a,P, 1 ) = 



4tt 
21 + 1 



(6.57) 



we finally obtain that the k coefficients are independent 
of the index m: 



K i (Al) = 0',M;^o|j,/u). 



(6.58) 



For (J, = j, on e find s nf ^ = t[ /z '. Indeed, according 



.1/2 



to the definition ( 5.15| ) of the SU(2) coherent states, the 
function Q(8, <fi) = (j, n\p\j, n) coincides with the proba- 
bility Pj(0, 4>) to find the displaced system in the highest 
spin state \j,j). It is not difficult to see that the proba- 
bility p_j(Q, (j)) to find the displaced system in the lowest 
spin state \j, — j) is e qual to Q(8 + n, (j)). As an applica- 
tion of the relation ( 6.46| ), we also obtain the following 
expression for the SU(2) Wigncr function in terms of the 
measured probability p M (n), 



2j 



1=0 



(4tt)- 1 (2Z + 1) 
(j,p;l,0\j,p) j v 



/ dn'Pi(n-n') Pfi (n'), 

(6.59) 



where Pi(x) are the Legendre polynomials. 



D. Informational completeness and unsharp 
measurements 



When the question of the state reconstruction arises, 
it is understood that the set of measurements one makes 
on an ensemble of identically prepared systems should 
give complete information about the quantum state. In 
particular, if one measures expectation values of some ob- 
servables, it is natural to ask how many such observables 
are needed to characterize completely the state of the 
system. In this sense a set of observables, whose expec- 
tation values are sufficient to reconstruct the quantum 
state (or, equivalently, to distinguish between different 
states), can be considered as informationally complete. 
A formal definition is as follows (6|] : A set of bounded 
operators 21 = {A} on H is said to be informationally 
complete if for density matrices p, p' the equality of ex- 
pectation values, 



Tr (pA) = Tr (p'A) MA e 21, 



(6.60) 



implies p = p' . 

The informational completeness of positive operator- 
valued measures covariant with respect to Heisenberg- 
Weyl, affine, and Galilei groups was recently discussed in 
Ref. [g3|. This subject was shown (6^] to be of importance 
not only in quantum mechanics but also in signal pro- 
cessing where a problem exists of extracting information 
from nonstationary signals and images. Another inter- 
esting feature is that both in quantum mechanics and in 
signal processing the phase-space formulation is of great 
importance for approaching this kind of problems. 
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It would be interesting to analyze the results of the 
present paper from the point of view of informational 
completeness. First, it is evident from the expansion 

P = Y,^{PDI)D V (6.61) 

that the orthonormal set of the tensor operators 

is informationally complete. Correspondingly, the set 
{A(f2; s) | Qe X} of the SW kernels for each s is also 
informationally complete. This fact is reflected by the 
inverse Weyl rule written as 

p= I d/z(fi)Tr[pA(fi;s)]A(fi;-s). (6.62) 
Jx 

In other words, the density matrix can be uniquely re- 
constructed from its s-parameterized QPD P(Q; s) — 
Tr [pA(f2; s)]. From the practical point of view, the im- 
portant thing is the informational completeness of the 
set {r u (f2) | fi G X} of the displaced projectors for 
any \u) G H. This fact was formally proved in Ref. 
p3| . Here, we presented a simple algorithm (based on 
the method of harmonic expansion) for the reconstruc- 
tion of the density matrix from the measurable probabil- 
ities p u (Q) — Tr [pr„(f2)]. This reconstruction procedure 
clearly implies the informational completeness of the set 

{r u (fi) \fiex}. 

One of the useful features of the method of displaced 
projectors is the ability to take into account in a sim- 
ple way the unsharpncss of a realistic measurement. Of 
course, it is impossible in practice to make a completely 
accurate displacement to a specified point f2 on the phase 
space. For example, one should take into account the 
phase and intensity fluctuations of a classical microwave 
source that displaces the quantum state of the radiation 
field in a cavity, or instabilities of a classical driving field 
that displaces the motional state of a trapped ion. Simi- 
larly, the so-called coarse graining problem arises in radar 
analysis due to frequency instabilities of the test signal 
or uncertainties in timing of signal initiation. As a re- 
sult, the probabilities should be integrated over 
the variation range. This yields the expectation value of 
the localization operator defined by ]63|| 

Z u (f)= [ dn(Sl) f(Sl)T u (Sl), (6.63) 
Jx 

where /(f2) is a localization function. Z u (f) has a purely 
discrete spectrum, is bounded when / G L k (X, n), k > 1, 
and is self-adjoint when / is real. In particular, let B be 
a region (more specifically, a Borel set) in X, with the 
characteristic function xs(^) that equals 1 for fi G B 
and otherwise. Taking /(£!) = xs(Sl), one obtains 

z u (B) = [ d^(n)T u (ri). (6.64) 

JB 

It is not difficult to see that the informational complete- 
ness of the set {r„(fi) | £1 e X} of displaced projec- 
tors implies the informational completeness of the set 



{Z U (B) | B G Borel sets of X} of localization operators. 
Therefore, in the case of realistic unsharp measurements, 
the localization operators may be conveniently used for 
analysis and reconstruction of quantum states or elec- 
tronic signals and images. 



VII. CONCLUSIONS 

In the present paper we propose a simple algorithm for 
constructing the SW kernels which implement the lin- 
ear bijective mapping between Hilbert-space operators 
and phase-space functions for physical systems possess- 
ing Lie-group symmetries. The constructed kernels are 
manifestly covariant under the action of the correspond- 
ing dynamical symmetry group and satisfy the traciality 
condition which ensures that quantum expectations can 
be represented by statistical-like averages over the phase 
space. Adding the noncommutative twisted product that 
equips phase-space functions with the algebraic structure 
of quantum operators, an autonomous phase-space for- 
mulation of quantum mechanics is developed. 

It turns out that the concept of phase space naturally 
emerges in the majority of schemes proposed for the re- 
construction of quantum states as well as in the standard 
methods of signal analysis. In particular, we focus on the 
method based on measurements of displaced projectors 
and develop its general group-theoretic description. We 
do so by applying the same technique of harmonic ex- 
pansions on the phase space that was used for the con- 
struction of the SW kernels. The problem of the state 
reconstruction is also approached using the concept of in- 
formational completeness, and the role of localization op- 
erators in describing realistic measurements is discussed. 
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